Counterion correlations and attraction between like-charged macromolecules 
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A simple model is presented for the appearance of attraction between two like charged polyions 
inside a polyelectrolyte solution. The polyions are modeled as rigid cylinders in a continuum dielec- 
tric solvent. The strong electrostatic interaction between the polyions and the counterions results 
in counterion condensation. If the two polyions are sufficiently close to each other their layers of 
condensed counterions can become correlated resulting in attraction between the macromolecules. 
To explore the counterion induced attraction we calculate the correlation functions for the con- 
densed counterions. It is found that the correlations are of very short range. For the parameters 
specific to the double stranded DNA, the correlations and the attraction appear only when the 
surface-to-surface separation is less than 7 A. 

PACS number(s); 61.20.Qg, 87.15.Aa, 87.15.-v 



I. INTRODUCTION 

In the last few years a new phenomenon has at- 
tracted attention of the community of soft condensed 
matter physicists — appearance of attraction between 
like charged macromolecules in solutions containing mul- 
tivalent ions. The problem is particularly fascinating be- 
cause it contradicts our well established intuition that 
like charged entities should repel The fundamen- 

tal point, however, is that the electrolyte solutions are 
intrinsically complex systems for which many body in- 
teractions play a fundamental role. 

The attraction between like charged macromolecules 
is important for many biological systems. One particu- 
larly striking example is provided by the condensation of 
DNA by multivalent ions such as Mn^"*", Cd^"*" and vari- 
ous polyamines . This condensation provides an an- 
swer to the long standing puzzle of how a highly charged 
macromolecule, such as the DNA, can be confined to a 
small volume of viral head or nuclear zone in procaryotic 
cell. Evidently, the multivalent ions serve as a glue which 
keeps the otherwise repelling like-charged monomers in 
close proximity Q. In eukaryotic cells, the cytosol is tra- 
versed by a network of microtubules and microfilaments 
— rigid chains of highly charged protein (F-actin) — 
which in spite of large negative charge agglomerate to 
form filaments of cytoskeleton The actin fibers are 
also an important part of the muscle tissue, providing a 
rail track for the motion of molecular motor myosin. 

Although the nature of attraction between like charged 
macromolecules is still not fully understood, it seems 
clear that the attractive force is mediated by the multiva- 
lent counterions ||9|-p^ . A strong electrostatic attraction 
between the polyions and the oppositely charged multiva- 



lent counterions produces a sheath of counterions around 
each macromolecule. The condensed counterions can be- 
come highly correlated resulting in an overall attraction. 
It is important to note that the complex formed by a 
polyion and its associated counterions does not need to 
be neutral for the attraction to arise. Under some con- 
ditions the correlation induced attraction can overcome 
the monopolar repulsion coming from the net charge of 
the complexes. 

Recently a simple model was presented to account for 
the attraction between two lines of charges ]T6|"p^ . Each 
line had Z discrete uniformly spaced monomers of charge 
— g, and n < Z/a condensed counterions of charge aq free 
to move along the rod. The net charge of such a polyion- 
counterion complex is Q,,{{ = —{Z — an)q < 0. Neverthe- 
less, it was found that if n > Z/2a and a > 2, at suffi- 
ciently short distances, the two like-charged rods would 
attract It was argued that the attraction resulted 

from the correlations between the condensed counterions 
and reached maximum at zero temperature. If n < Z/2a 
the force was always found to be repulsive. 

Clearly, a one dimensional line of charge is a dramatic 
oversimplification of the physical reality. If we are inter- 
ested in studying the correlation induced forces between 
real macromolecules their finite radius must be taken into 
account Thus, a much more realistic model of 

a polyion is a cylinder with a uniformly charged back- 
bone 1^,^ or with an intrinsic charge pattern [pT|, 22 
as, e.g., the helix structure of DNA molecule. Further- 
more, the condensed counterions do not move along the 
line, but on the surface of the cylinder. Unfortunately, 
these extended models are much harder to study analyt- 
ically. 

In this paper we explore the effects of finite polyion di- 
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ameter on the electrostatic interactions between the two 
polyions using Monte Carlo simulations. We find that 
the finite diameter and the associated angular degrees 
of freedom of condensed counterions significantly modify 
the nature of attraction. Thus, although there is still a 
minimum charge which must be neutralized by the coun- 
terions in order for the attraction to appear, this fraction 
is no longer equal to 50% as was the case for the line 
of charge model. We find that the critical fraction de- 
pends on the valence of counterions and is less than 50% 
for a > 2. For monovalent counterions no attraction is 
found. The crystalline structure of the condensed coun- 
terions, as first suggested by simulations of Gronbech- 
Jensen et al. ^ and Refs. [^2|jl|,|l|,|§ , is also not very 
obvious. In particular we find very similar distributions 
of condensed counterions in the regime of attractive and 
repulsive interactions. 

The structure of this paper is as follows. The model 
and the method of calculation are described in section ||. 
In section III , we present the results of the simulations. 



The conclusions are summarized in section [V 




FIG. 1. Schematic depiction of the model: two polyions of 
Z negative charges (small solid circles at the centers of rings) , 
with radius R and n condensed counterions (large solid cir- 
cles) each, are separated by distance d. The counterions are 
free to move along the Zm positions (open circles) fixed on 
the rings located around each of the Z monomers. In this 
sketch, Z = A, n = 2 and m = 4. 



monomers, each polyion has n < Z/a condensed coun- 
terions with valence a and charge aq, which are con- 
strained to move on the surface of the cylinder. To lo- 
cate a condensed counterion it is necessary to provide its 
longitudinal position, z (0 < z < Z), and the transver- 
sal angle, {0 < 9 < 2tt). To simplify the calculations, 
the angular and the longitudinal degrees of freedom are 
discretized, see Fig. |l|. The surface of the cylinder is sub- 
divided into Z parallel rings with a charged monomer at 
the center of each ring. Each ring has m sites available to 
the condensed counterions, see Figs. |^ and ^. The hard- 
core repulsion between the particles requires that a site 
is occupied by at most one condensed counterion. The 
two polyions are parallel, with the intermolecular space 
treated as a uniform medium of dielectric constant e. 




FIG. 2. The cross sectional view of the adjacent rings 
on the two polyions. The sites are labeled with integers 
J = 0, 1, . . . , m — 1 in the counter-clockwise direction. The 
angle between two consecutive sites is A9 = 2n/m. Here we 
show the angles Ol of site z = 1 on the polyion 1, and 9^ of 
site i = 3 on the polyion 2. 

We introduce occupation variables n*^ for the two 
polyions, so that i = 1,2,..., Z{m + 1) and k = 1,2. 
Thus, = 1 if the I'th site of the k'th polyion is occu- 
pied by a particle of valence ai = { — 1, a} (negative core 
charge or counterion of valence a, respectively), other- 
wise = 0. Note that the core charge is always "oc- 
cupied" , while the counterions are free to move between 
the Zm ring sites of each polyion. The Hamiltonian for 
the interaction between the two polyions is. 



(1) 
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II. MODEL AND METHOD 

A. The model 

The DNA model considered here is an extension of the 
one proposed earlier by Arenzon, Stilck and Levin ||l6|,|l^ . 
A similar model has been recently discussed by Solis and 
Olvera de la Cruz |Q . The polyions are treated as par- 
allel rigid cylinders of radius R and Z ionized groups, 
each of charge —q, uniformly spaced — with separation 
b — along the principle axis. Fig. |l|. Besides the fixed 



with the i ^ j ■ All lengths are measured in units of h (for 
DNA, b — I.7A). The dimensionless quantity, ^ = I3q'^ /be, 
is the Manning parameter, which for DNA is ^ = 4.17. 
The partition function is obtained by tracing over all the 
possible values of n\ consistent with the constraint of 
fixed number of counterions per polyion, 



(2) 



Clearly, this is a very crude model of the interaction be- 
tween two macromolecules in a polyelectrolyte solution. 
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The molecular nature of the solvent is ignored. Also the 
number of condensed counterions is fixed instead of being 
dependent on the separation between the particles. Nev- 
ertheless, we believe that this simple model can provide 
some useful insights for the mechanism of attraction in 
real polyelectrolyte solutions. 
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FIG. 3. Mean site occupation for the adjacent rings on the 
two polyions located at z = 9. Each polyion has Z — 20 and 
n = 7. The distance between the centers of rings is d/b — 32.8 
(squares) and d/b = 16.65 (circles). The solid line corresponds 
to the sites on the first polyion and the dashed lines are for the 
second polyion. Note the almost perfect symmetry between 
the two macromolecules. 



B. The observables 

We are interested in statistical averages of observ- 
ables such as the energy and the force between the two 
polyions. Furthermore, to understand the nature of the 
interaction between the two macromolecules it is essential 
to study the correlations between the condensed counte- 
rions on the two polyions. 

The force is obtained from the partition function 
Eq. (H), b/3F — — V(lnQ). From symmetry, only the 
y-component is different from zero. 

For finite macromolecules the symmetry between the 
two polyions can not be broken p8| . Hence it is impos- 
sible to produce a true crystalline order in a finite sys- 
tem at non-zero temperature. Since within our simplified 
model the two polyions have exactly the same number of 
condensed counterions, the average angular counterion 
distribution {n^{z,6^)) must be symmetric with respect 
to the mid-plane y = d/2, see Fig. |l|. The angle 9'^ labels 
the site i on polyion k, see Fig. ||. Thus, n^{z, 6*1) denotes 
the occupation variable for the site 3 on the ring z located 
on polyion 2, with an angle of 3(27r/rn). Indeed, Fig. ^ 
shows that the density profiles are completely symmetric 
(up to fiuctuations). In spite of this symmetry it is pos- 
sible for the counterions on the two polyions to become 
highly correlated. Clearly, the strength of these correla- 
tions will depend on the product ^a^ and the separation 



between the two macromolecules. Considering Fig. g, it 
is evident that if the site 2 on the first polyion is oc- 
cupied, the likelihood of occupation of the site 8 on the 
second polyion will be reduced. 

To explore the nature of electrostatic correlations, we 
define a counterion-hole correlation function between the 
adjacent rings on the two polyions. 



cll = {nl{z,9l)[l-n'^{z,9])]) 

-{nliz,9l)){[l-n]iz,9])]) 



(3) 



Here (• • •) denotes the ensemble average. This function 
should be non-zero when sites on the two polyions are 
correlated, that is if one is occupied by a condensed coun- 
terion there is an increased probability of the second be- 
ing empty. 



C. Simulations 



To calculate the force between the two polyions, we 
have performed a standard Monte-Carlo (MC) simula- 
tion with the usual Metropolis algorithm p^]. First, one 
counterion on polyion 1 is randomly chosen and displaced 
to a vacant position on the same polyion. This move is 
accepted or rejected according to the standard detailed 
balance criterion |2^ . We do not permit exchange of par- 
ticles between the polyions. Next, the same is done for 
polyion 2. In one Monte-Carlo step (MCS) all 2n con- 
densed counterions on the two polyions are permitted to 
attempt a move. 

The long-ranged nature of the Coulomb interaction re- 
quires evaluation of all the pair interactions in Eq. (|l|) 
at every MCS. Due to the limited computational power 
available to us, we have confined our attention to rela- 
tively small systems with Z < 100 and m = 10. We have 
checked, however, that for m — 10 the force has already 
reached the continuum limit and did not vary further 
with increase of m. Also we note that the "thermody- 
namic limit" is reached reasonably quickly, so that there 
is a good collapse of data already for Z > 50, see Fig. 
2000 MCS served to equilibrate the system, after which 
500 samples were used to calculate the basic observables, 
namely, the mean force and energy. To obtain the corre- 
lation functions, 5000 samples were used with 5000 MCS 
for equilibration. 



III. RESULTS AND DISCUSSION 



The simulations were performed for ^ = 4.17 and 
R/b = 8.2, relevant for DNA. For monovalent counteri- 
ons the simulation results indicate that the force is purely 
repulsive. This is in complete agreement with the ex- 
periments , which do not find any indication of DNA 
condensation for monovalent counterions. 
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FIG. 4. Mean force in j/-direction versus distance d/b be- 
tween the two DNA molecules R/b = 8.2, Z ^ 20 and 
^ — 4.17. The symbols indicate the number n of divalent 
{a = 2) condensed counterions. 

For divalent counterions the force between the two 
complexes can becomes negative, indicating appearance 
of an effective attraction. Fig. ^ The range of attraction 
is larger than was found for the one dimensional line of 
charge model, Ref. 
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FIG. 5. The surface-to-surface separation below which the 
force becomes attractive, as a function of the number n of 
condensed counterions, for valences a — 2 (full) and 3 (open) , 
^ = 2.283 (circles) and 4.17 (squares), with Z ranging from 
50 to 100. 

Within the Manning theory H 88% of the DNA's 
charge is neutralized by the divalent counterions. How- 
ever, there are indications that even a larger fraction 
of DNA's charge can become neutralized by the multi- 
valent ions if the counterion correlations are taken into 
account [Q. In this case the interaction is purely at- 
tractive, with the rang e of about d/b « 20 or 34 A (7 A 
surface-to-surface). Fig. ^. 
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FIG. 6. Snapshots of two equilibrium configurations for (a) 
d/b — 32.8 and (b) d/b — 16.8, for two polyions with Z — 60 
and n — 18. 

A minimum number of condensed counterions is nec- 
essary for attraction to appear. In Fig. ^ we present 
the surface-to-surface separation, (do — 2i?), below which 
the force between the two complexes becomes negative 
(attractive), as a function of the number of multivalent 
counterions. For the case of DNA with divalent counte- 
rions a = 2, the attraction appears only if 40% of the 
core charge is neutralized. For a = 3 this fraction de- 
crease to 30%. Furthermore, decrease in the value of the 
Manning parameter, ^, increases the minimum number 
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of condensed counterions necessary for the attraction to 
appear. This is fully consistent with the fact that the 
attraction is mediated by the correlations between the 
condensed counterions. Since raise in temperature tends 
to disorganize the system, the state of highest correlation 
between the condensed counterions corresponds to T = 
or ^ = oo. 

The surface-to-surface distance at which the attrac- 
tion first appears tends to zero as the number of con- 
densed counterions is diminished. We find do — 2R ^ 
{fj, — iJLc{oi)Y , where the average counterion concentra- 
tion is /i = n/Z and the critical fraction /Zc depends on 
the valence of condensed counterions a. From Fig. || it is 
evident that v = 1. This should be contrasted with the 
line of charge model Ref. for which v = 1/3. 
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FIG. 7. Site-site correlation functions, Eq. for the i'th 
site on the ring z; = 9 of the first polyion, with the j'th site 
on the z = 9 ring of the second polyion. We consider only 
the correlations of the four "inner sites" of the first polyion 
i = 1,2, 3, 4, (the four curves), with all the sites of the 2 = 9 
ring of the second polyion, j — 1...9. The total number of 
sites per ring is m = 10. 9^ indicates the angular position 
of the j'th site (in units of 2n/m) on the ring of the second 
polyion. The parameters are as in Fig. ^ with n — 7 and the 
separation between the two macromolecules is d/h = 32.8. 
The graph shows that at this distance there are almost no 
correlations between the condensed counterions. From Fig. |^ 
we also see that the effective force is repulsive. 

In Fig. H we show two snapshots of the characteris- 
tic equilibrium configurations for (a) d/b = 32.8 and (b) 
d/b = 16.8. Looking at this figures it is difficult to see 
something that would distinguishes between them, both 
appear about the same. There is no obvious crystal- 
lization or transversal polarization suggested in previous 
studies 20|. Yet, the case (a) corresponds to the re- 
pulsive, while the case (b) corresponds to the attractive 
interaction between the polyions. To further explore this 
point, in Figs. |^ and || we present the site-site correla- 
tion function, Eq. (H), for macromolecules with Z = 20 



and n = 7. For d/b = 32.8 the surface-to-surface dis- 
tance between the two polyions is sufficiently large for 
their condensed counterions to be practically uncorre- 
lated. Fig. |^. On the other hand, for d/b ~ 16.65 strong 
correlations between the condensed counterions are ev- 
ident. Fig. ||. The Fig. || shows that the sites two and 
three on the first polyion are strongly correlated with the 
sites seven and eight on the second polyion, respectively. 
It is these correlations between the adjacent sites on the 
two polyions which are responsible to the appearance of 
attraction between the two macromolecules when they 
are approximated. Fig. ^. 
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FIG. 8. Site-site correlation functions for separation 
d/b — 16.65, where, according to Fig. ^, there is attraction. 
Note that sites 2 and 3 (open diamond and triangle, respec- 
tively) on polyion 1 are strongly correlated with sites 8 and 7 
on polyion 2, respectively 



IV. SUMMARY 

We have presented a simple model for polyion-polyion 
attraction inside a polyelectrolyte solution. It is clear 
from our calculations that the attraction results from 
the correlations between the condensed counterions and 
reaches maximum for T — 0. The thermal fluctuations 
tend to diminish the correlations, decreasing the ampli- 
tude of the attractive force. Consistent with the exper- 
imental evidence, the attraction exists only in the pres- 
ence of multivalent counterions. Our simulations demon- 
strate that a critical number of condensed counterions is 
necessary for the appearance of attraction. The fraction 
of bare charge that must be neutralized for the attrac- 
tion to arise depends on the valence of counterions. The 
larger the valence, the smaller the fraction of the bare 
polyion charger that must be neutralized for the attrac- 
tion to appear. This result should be contrasted with 
the line of charge model ||l^ for which the critical frac- 
tion was found to be equal to 50%, independent of the 
counterion charge. 
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